Locality in Quantum and Markov Dynamics on Lattices and Networks 
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We consider gapped systems governed by either quantum or Markov dynamics, with the low-lying 
states below the gap being approximately degenerate. For a broad class of dynamics, we prove that 
ground or stationary state correlation functions can be written as a piece decaying exponentially in 
space plus a term set by matrix elements between the low-lying states. The key to the proof is a 
local approximation to the negative energy, or annihilation, part of an operator in a gapped system. 
Applications to numerical simulation of quantum systems and to networks are discussed. 



PACS: 05.30.-d,03.67.Lx 

The relation between energy and length scales is cen- 
tral to our understanding of physics. We intuitively asso- 
ciate low energies with long wavelengths. For many-body 
systems, at a quantum critical point, where the energy 
gap vanishes, we expect to see long-range correlations 
[1]. Conversely, experience teaches us that a gap in a 
quantum system implies a finite correlation length. 

Such a result is well known for non-interacting systems; 
for example, a defect in a perfect crystal may give rise to 
localized modes lying within a band gap. The basic result 
shown in [2] is that for a large class of short-range many- 
body Hamiltonians, this result still holds: if there is a 
gap between the ground and first excited states, then the 
connected correlation function of two operators decays 
exponentially in space with correlation length £ bounded 
by a characteristic velocity divided by the gap, AE. We 
note that this contrasts with the possibility of having an 
infinite entanglement length in a gapped system, as is 
studied in the field of quantum computation [3]. 

In this Letter, we extend this result to a more gen- 
eral class of systems. We consider the case of a quantum 
system with some number of almost degenerate low en- 
ergy states, all within energy AEi ow of the ground state 
energy, with the rest of the spectrum having an energy 
at least AE above the ground state. Then, we prove 
that the correlation functions include an exponentially 
decaying piece with correlation length £, plus a piece 
which involves matrix elements between the states be- 
low the low-lying states. The results in this Letter apply 
to finite range lattice Hamiltonians, with some techni- 
cal conditions required to bound the group velocity on 
the lattice, as discussed more below. This set of systems 
includes short-range lattice spin systems, lattice fermion 
systems, and lattice hard-core boson systems. This result 
provides a general proof of Kohn's idea of "nearsighted- 
ness" [4] for this class of systems, and thus may have 
important applications in O(N) methods for simulating 
quantum systems. 

Proofs of some of the results used in this Letter can 
be found in [2] , and additional proofs will be given in [5] . 
After giving the basic results, we discuss a wide variety 
of applications: to systems with a band structure; to 



classical Markov processes; and to systems on general 
networks. 

Consider a given connected correlation function: 
(AB), where A,B are operators and (...) denotes the 
ground state expectation value. Using a spectral rep- 
resentation, we have (AB) = Y1e<ae Aw-Bio + 
J2e >AE AoiBio, where i represents different intermedi- 
ate states and is the ground state. Define (Ai ow )ij = 
Ay if both Ei and Ej are less than or equal to AEi ow , 
while (Ai ow )ij — otherwise. 

In this Letter we determine which correlation functions 
may be long-range, given the structure of low energy 
states. The basic result is that for a pair of operators 
A, B separated by distance I, for vanishing AEi ow , 



(AB) 



{A low ,B low }) +C?(exp[-Z/£] 



(1) 



A quantum Ising system with a transverse field is 
a good example system to apply this result: H = 
Jj2d{i,j)<i SfSj + B J2i "Sf> where Sf are the spin op- 
erators on site i for a = x,y, z and where d(i,j) is some 
metric on the lattice. In the paramagnetic phase (B/J 
sufficiently large), the system has a unique ground state 
with a gap to the rest of the spectrum and so Eq. (1) 
implies that all connected correlation functions decay ex- 
ponentially. In the ferromagnetic phase (B/J sufficiently 
small) , the system has two almost degenerate low energy 
states, and again has a gap to the rest of the spectrum. 
Operators such as Sf have long-range correlations in the 
ferromagnetic phase, due to matrix elements of these op- 
erators between the two low-lying states. Correlation 
functions of operators which do not couple the low-lying 
states, such as energy-energy correlation functions, are 
exponentially decaying. 

States Above the Gap — Define 

1 f°° 1 
A ± = 7 F dtA(t)——, (2) 

where A(t) = ex.p(iHt)Aexp(-iHt). Then, (A+) tJ = 
AijO(Ei — Ej), where G is a step function: Q(x) = 1 for 
x > and S(x) = for x < 0, while 6(0) = 1/2. Thus, 
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A~ includes only the negative energy (positive frequency) 
matrix elements of A. We may define a similar for 
any operator 0. 

Define A high = A - A iow Then, 



(AB) = (A low B low ) + (A- tgh B) 
= (A low Bi ow ) + {[A Mgh ,B]), 



(3) 



since {BA^ igh ) — 0. Here, we may view A + as a creation 
operator and A~ as an annihilation operator. 

The basic idea of the proof in this Letter is that we 
will find an approximation to A~ such that (1): the er- 
ror involved in the approximation is small and (2): the 
commutator of the approximate operator with B is small. 
The approximation is defined by: 



A(t) = A(t)exp{-(tAE) 2 /(2q)}, 



2tt/ 



dtA(t) 



1 



(4) 
(5) 



where q will be chosen later. In general, we may define 
0,0 + for any operator O. 

Then {A + )ij = A i] Q q (E l - Ej), where 



Jo 



do/ 
~2tT 



{^2Trq/AE) exp[-g(w - J) 1 /{2AE 2 )]. 

(6) 



The function 9 9 is equal to an error function; for large q 
it approximates the step function 0. 

Finite Velocity — Since A and B are separated in space, 
they commute. Suppose operators A, B are separated in 
space by a distance I, meaning that I is the shortest dis- 
tance between any two sites i, j, such that some operator 
on site i appears in A and some operator on site j ap- 
pears in B. Then, one expects that the commutator of 
the operators A(t) and B will be small for some range 
of times t, with \t\ less than or equal to I divided by a 
characteristic velocity of the system. 

Specifically, we consider any lattice Hamiltonian 7i 
which can be written as a sum H = ^ i H l , where i ranges 
over lattice sites and where we require that: (1) the com- 
mutator [H % , O] = for any operator O which acts only 
on sites j with d(i,j) > R, where R is the range of the 
Hamiltonian; and (2) the operator norm \\Tt l \\ < J, for 
all i, for some constant J. Then, it was shown in [6,2] 
that there exists a function g(t,l), which depends on J, 
R, and the lattice structure, such that 



IP(t),fl(0)]||< 



B||J>(Mi). 



(7) 



where the sum ranges over sites j which appear in opera- 
tor B, and lj — d(j, i) is the distance from j to the closest 
site i in the operator A. It was shown that there exists 
some constant c\ such that g(c\l,l) is exponentially de- 
caying in I for large I [2]. Further g is symmetric in t so 



that g(t,l) = g(-t,l) and g(t,l) < {t/f)g(f,l) for t < t' 
and t, t' > 0. Thus, g{t,l) is monotonically increasing in 
t for t > 0. Hence, for t < Cil, \ \[A(i), B(0)]\ \ is exponen- 
tially decaying in I for large I. We define v = c^ 1 as a 
velocity of the system. We will define 



C(t) = 5>(Mi), 



(8) 



where the sum ranges over sites j which as above. In 
a later section, we will discuss the applicability of this 
result to systems on general graphs, or networks. 

As a side note, this bound on the operator norm of 
the commutator may be viewed as a non-relativistic ana- 
logue of the concept of micro-causality in relativistic field 
theories [7]. In relativistic theories, operators commute 
outside the light-cone, so that if c is the speed of light, 
then [A(t), -B(O)] = for t < l/c. In a non-relativistic 
theory, the commutators do not exactly vanish, but we 
are able to find a velocity v such that the commutator is 
exponentially small for t < I /v. 

Commutator — To evaluate the commutator in Eq. (3), 
we use 

{[A Mgh ,B]) = ([A-,B]) + {{[A- high ,B]) {[A- high ,B])) (9) 

-<fc,.*]>. 

Here, A^ igh is defined by starting with A^igh, and then 
multiplying by exp[— (tA) 2 /(2q)], following Eq. (4), and 
finally taking the negative energy part. Thus, A^ igh + 

Aow = A ■ 

First we bound the first term on the right-hand side of 
Eq. (9). We have 



\([A-,B})\ (10) 
1 





/ 

J\t 








j 

J\t 


< 


1 

27 



\t\<d 



\t\>dl 



At C M-{t^Ef/{2q)]{[A{t) 1 B]) — 
dtcM-(t^E) 2 /(2q)]([A(t),B]) 



it + e 
1 



it + e' 



AE Cl l 



)■ 



In Eq. (10), for \t\ < c\l, we have 
used eM-(t^E) 2 /(M([Mt),B)]) < \\[A(t),B)]\\ and 
Eq. (7), while for 

|i| > dl we have used exp[—(tAE) 2 /(2q)]([A(t), B}]) < 
2exp[-(tAE) 2 /(2q)]\\A\\\\B\\, and then we have per- 
formed the integrals using these bounds. 

Next, we consider the second pair of terms on the right- 
hand of Eq. (9). From Eq. (6) it follows that, for \u>\ > 
AE, \e q (uj) - e(uj)\ < exp[-g/2]/V2lFg. Thus, if * ) is 
the ground state wavefunction, \{9 A^ igh - {^oA Mgh \ < 
\\A\\exp[-q/2]/y/TFq, and also \A^ igh * ) - A^ flh * >| < 
\\A\\exp[-q/2]/ y / ! Znq. Then, 



2 



\([An igh ,B})-([A- igh ,B})\<2\\A\\\\B\ 



exp[-q/2] 



(11) 



Physically, Eq. (11) is a kind of uncertainty relation: the 
time integral in Eq. (4) extends over a time of order 
y/qAE^ 1 , and thus provides an approximation to ener- 
gies of order AE/^/q. 

Almost Degenerate Low Energy States — 
Suppose AEi ow is very small compared to AE. 

< AEi ow , Qq(cj) is close to 1/2: 



Thus, for AE low = 0, we find that \([Ai ow , Bi ow ])\ is ex- 
ponentially decaying in I. 

Operators at Different Times — We finally extend the 
result Eq. (13) to correlation functions (A(— ir)B(0)), 
with r real and r > 0. Define 



A ± {±ir) 



2tt/ 



dti(t) 



±it + T 



(15) 



Then, for ui 

|0,(w) - l/2\ < Jq~/2^(w/AE) 



\{[Aiow,B])\ < 2\\Ai ow \\\\Bi ow 
from Eqs. (3,9,10,11) above, it follows that 



Thus, \([A- ow ,B}) - 
\\y/qj2^(u/AE). Then, 



< 2 



B\ 



MB) 
'C(ci0 



{Ai ow Bu 
2tt ' 2-kAEciI 



+ 



)--([A low ,B low })\ (12) 
exp[~( Cl lAE) 2 /(2q)} 



by —it + t. In this case, we find that Eq. (10) 
still holds as a bound for \([A~ (—ir), B})\. One may 
also show that, for r < q/AE, \([A^ igh (—iT), B]) - 
C M+r 2 ^/{2q)]{[A Mgh (-ir) 1 B])\ < 2||A||||B|| e ^/ 2 . 

With the given q = c\lAE the above bounds show that 
for r < ciZ, \(A(-ir)B(0)) - (A low {-ir)B low {0))\ < 

e -^AS/(2 Ca 0[_L 2 || A |||| jB || C - (ci/) + (16) 

2tt 



+ 



cxp 



+ 2\\A low \\\\B low \\^/q~/2^(AE low /AE). 



Picking q = cilAE, we have 



\(AB)-(-{A, 



< 




(13) 



B||[2C( Cl 0/(27r) 

+2| \A low 1 1 1 \B low 1 1 v / c 1 lAE/2ir(AEi ow /AE) . 

The first term on the right-hand side of Eq. (13) is expo- 
nentially decaying in I with some correlation length ^c- 
The terms in parenthesis are exponentially decaying with 
correlation length 2/(ciAE). Thus, we define the corre- 
lation length £ to be the minimum of £c and 2/(c\AE). 
For AEi ow taken to be zero, we have Eq. (1). 

Low Energy Commutator — We can also bound 
the expectation value ([Ai OW7 Bi ow ]). Define A — 
AE/y/2-Kq dt A{t). Then, we can bound the com- 
mutator: \([A°,B])\ < A£/v^Fg(jj t|<ciJ dt C(t) + 
\J w>Cll dt{[A(t),B]\) < \\A\\\\B\\^c 1 lAE/2nC(c 1 l) + 

2e - Cl lAE/2y 

Next, consider the difference ([A , B]) — ([Ai ow , Bi ow \) . 
This equals (([A^B) - ([A low , B low })) + ([A° high ,B]). 
Explicit computa- 
tion with a spectral representation gives \([Af ow ,B) — 
([A low ,B low })\ < 2\\A low \\\\B low \\(e-i( AE '°-/ AE ) 2 / 2 - 1), 
while \{[A° high ,B])\ < 2\\A\\\\B\\e- CllAE / 2 . Combining 

these bounds with the bound on |(L4°,.B])| gives 



I ([A 



low ■> ^low 



< ^lAE/2n\\A\\\\B\\C( Cl l) (14) 
+4||A||||B||e- Cl ' A£ / 2 
+2\\A low \\\\B low \\(e-^ AE -^ A ^ 2 / 2 1). 



\B\\e 



ctAEl/2 



\([A low (-iT),B])\. 



To bound the last term in the above Equa- 
tion, we use \\Ai ow (t) - Ai ow (0)\\ < tAEi ow \\Ai ow \\. 
Define z = (27T)- 1 / dt exp[-(tAE) 2 /(2q)}/(-it + 
t). Then, \{[Ai ow {-ir),B]) - z([A low , B low })\ < 
ir^Jdt tAE low | \A low 1 1 1 \\B low 1 1 exp[-(i AE) 2 / (2q)]/ {-it+ 
t) < ^/2qJ^\\A low \\\\B low \\(AE low /AE). Thus, for 
AE low = 0, \([Ai ow (-iT),B])\ is exponentially decaying 
in I, following Eq. (14). 

Band Structure — Another interesting application of 
the above techniques is to problems with a band struc- 
ture. Note that all of the results above can be general- 
ized to fermionic operators A, B by interchanging com- 
mutators and anti-commutators throughout. Suppose we 
have a free fermionic theory, with a spectrum which has 
two bands with a band gap of width 2AE. Then, we 
can shift the zero of energy so that the spectrum has 
some set of states with energy at most —AE and an- 
other with energy at least AE. Then, if A = ipj is the 
fermionic creation operator at some point i, we can de- 
fine an operator A~ which approximately projects ipj 
onto the lower band. At the same time, A~ will be expo- 
nentially localized around point i, so that ||{j4 _ ,0}|| = 
0(exp[-(c 1 lAE) 2 /(2q)])+0(C(c 1 l)), if the fermionic op- 
erator O acts only on a site j with d(x,j) = I. This pro- 
jection technique may be a useful way to compute the 
density matrix in these systems [8]: if the chemical po- 
tential is such that all states are filled up to zero energy, 
then p(i,j) = (4^) = ({A~,ipj}) « ({A-,^}). 

Markov Processes — The above locality results can also 
be carried over to systems which obey continuous time 
dynamics, following [2], where we have a transition ma- 
trix Tij and a probability pi of being in state i, so that 
dtPi = J2j TijPj . By conservation of total probability, 
we have z2 i T i j = 0, guaranteeing that T has at least 
one zero eigenvalue. Let the stationary state with zero 
eigenvalue have right eigenvector let Li be the vector 
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with Ii = 1 for all i; this is a left eigenvector of T with 
eigenvalue 0. 

Suppose all eigenvalues of T are real (this includes all 
systems for which the stationary state obeys detailed bal- 
ance). All eigenvalues of T are non-positive. Then, as- 
sume that there are some number of eigenvalues A, of 
T with > Xi > —Ai ow , while all other eigenvalues 
Xi have Xi < —A, with A > A; otu . For each quan- 
tity to be measured, A, B, define (A) — J2i^-iPi- 
We can introduce for each quantity a diagonal matrix 
given by An — A i} and Aij = for i ^ j. Then, 
(A(t)B(0)) = lUxp[-Tt]Aexp[Tt}Bp° = PA(t)Bp°. 
We have left off the indices on the vector I,p and on 
the matrices A, B, exp[±Ti]; the product is evaluated 
following the usual rules of matrix multiplication. We 
can continue these definitions of A(t) to imaginary time, 
and define A(it) = A{it) cxp[-(iA) 2 /(2<?)] while A+ = 
±$&tA(it){Tit + e)-\ 

Assume that T can be written as a sum of matrices 
T % , with finite interaction range R and bound ||T l || < J 
for all i. With these preliminaries, all of the above ma- 
nipulations can be carried out for Markov processes. We 
find, in particular, that (AB) = {{\/2){Ai ow , Bi ow }) + 
0(exp[-l/f}). 

Networks — Recently, systems on general graphs or 
networks have been much studied [9]. Example systems 
include the random graph and the small-world network 
[10]. Locality is an important question in these systems, 
for community detection, for example [11,12]. In a small- 
world, the presence of long-range jumps can completely 
destroy locality above some length scale so that all criti- 
cal behavior is mean-field [13]. 

The study of quantum systems on networks is a bit 
unphysical, but Markov processes are extremely natural 
to study on networks. A good example is the contact 
process for epidemic spreading [14]. 

Following [2], the bound Eq. (7) holds for any graph 
for which all sites have a bounded number of neighbors 
(counting a neighbor as any other site within distance 
R) and for which |T 4 | < J for all sites i. This includes 
a large class of interesting networks, such as the small- 
world network. However, another large class of interest- 
ing networks, the scale- free networks, have an unbounded 
coordination number. In many of these cases, Ti. (or T) 
can be written as a sum of operators lying on bonds which 
have bounded operator norm. Then, we can add to the 
graph a set of additional sites which lie on these bonds 
between sites on the original network, calling one of new 
these sites if it lies on the bond between i and j. 

Then, we can let H = , with \\H^\\ < J. 

However, we still must deal with the unbounded num- 
ber of neighbors. The range of one of these bond Hamil- 
tonians is such that (i,j) is within range R of (i,k) for 
all j, k which neighbor i, so that we cannot bound the 
number of neighbors within range of a given bond 



However, if the given scale-free networks is loopless, in- 
spection of the proof in [2] shows that Eq. (7) still holds. 
Thus, we claim that for networks with either bounded co- 
ordination number and bounded Ti. 1 or loopless networks 
with bounded H 1 ^ , the result Eq. (1) is valid. 

As an application, consider a contact process in the 
endemic phase on a small- world network. The Markov 
process governing the dynamics has a gapped transition 
matrix away from the critical point. Then, a correlation 
function of the number of infected individuals on a site 
i (this number is zero or one) with the number on a site 
j decays exponentially as exp[—d(i,j)/£\ where d(i,j) is 
the shortest path distance (including long-range links) 
from i to j. 

Discussion — The result we have shown is expected on 
physical grounds. In a gapped system, the correlation 
function can be written as (1): an exponentially decay- 
ing piece, plus a piece determined by matrix elements be- 
tween the low lying states. In the quantum Ising system 
mentioned in the introduction, the two low-lying states 
are symmetric and anti-symmetric combinations of spin- 
up and spin-down states. Acting on the system with a 
local operator, such as Sf, changes between these two 
states, and so one may have long-range spin-spin correla- 
tion functions; these correlation functions converge expo- 
nentially to some constant at large separations. In other 
cases, such as a spin liquid or fractional quantum Hall 
system with a topological degeneracy and a gap to the 
rest of the spectrum [15], the matrix elements between 
the low-lying states are small for all local operators, so 
all correlation functions decay exponentially. 
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